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(57) Abstract: A magnetic resonance imaging method is proposed wherein a magnetic resonance image is reconstructed from mag- 
netic resonance signals from respective signal channels. More specifically, individual signal channels relate to respective surface 
coils which are employed as receiver antennas for the magnetic resonance signals. The magnetic resonance signals are acquired 
with sub-sampling of the k-space. Resampling on a regular square grid is performed, thus enabling fast Fourier transformation in 
the reconstruction of the magnetic resonance image. Furthermore, the reconstruction is carried out on the basis of the spatial sensi- 
tivity profile of the receiver antennas, i.e. of the surface coils, so as to separate contributions from different spatial positions in the 
sub-sampled magnetic resonance signals. Preferably, a spiral-shaped acquisition trajectory is followed in the k-space. 
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Magnetic resonance imaging method with sub-sampling 



The invention relates to a magnetic resonance imaging method for forming a 
magnetic resonance image, wherein 

magnetic resonance signals are acquired by means of receiving antennas via a 
plurality of signal channels, 

5 - which individual receiving antennas have respective sensitivity profiles. 

The invention also relates to a magnetic resonance system. 
A magnetic resonance imaging method and a magnetic resonance system for 
carrying out such a magnetic resonance imaging method are known from the article "Coil 
Sensitivity Encoding for Fast MRI" by. K.P. Prussmann et al. in Proceedings ISMRM (1998), 
10 579. 

The known magnetic resonance imaging method is known by the acronym 
SENSE method. This known magnetic resonance imaging method utilizes receiving antennas 
in the form of receiving coils. This magnetic resonance imaging method utilizes sub- 
sampling of the acquired magnetic resonance signals so as to reduce the time required to scan 

1 5 the k-space at a sampling density in the k-space for the desired field-of-view and over a 
region in k-space which is large enough for the desired spatial resolution of the magnetic 
resonance image. Notably the respective lines in the k-space along which scanning is 
performed are situated apart further in the k-space than is necessary for the desired spatial 
resolution. In other words, "lines are skipped" in the k-space. As a result of such "skipping of 

20 lines in the k-space", less time is required for the acquisition of the magnetic resonance 

signals. Receiving coil images are reconstructed on the basis of the sub-sampled magnetic 
resonance signals from the individual receiving coils. Due to the sub-sampling, the actual 
field of view is reduced so that back-folding or aliasing artefacts occur in such receiving coil 
images. The magnetic resonance image is derived from the receiving coil images on the basis 

25 of the sensitivity profiles, the aliasing artefacts in the magnetic resonance image thus being 
substantially or even completely removed. This unaliasing operation enlarges the magnetic 
resonance image to the desired field of view. 

It has been found in radiological practice the time required for the acquisition 
of the magnetic resonance signals need still be reduced considerably further. It has been 
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found that a substantial reduction of the acquisition time of the magnetic resonance signals is 
necessary notably for magnetic resonance imaging methods for the imaging with a high 
spatial resolution of fast moving parts of the anatomy, for example the beating heart of a 
patient under stress, but also for MR angiography. 
5 It is an object of the invention to provide a magnetic resonance imaging 

method wherein the acquisition time of the magnetic resonance signals is significantly shorter 
than the acquisition time required when the known SENSE technique is used. 

This object is achieved by means of the magnetic resonance imaging method 
according to the invention wherein the noise correlation between individual signal channels is 
10 represented by a noise correlation matrix, where 

the magnetic resonance signals are acquired with sub-sampling, 

regularly resampled magnetic resonance signals are resampled on a regular 
sampling grid from the acquired magnetic resonance signals, 

the noise correlation matrix is approximated by a block diagonal matrix or a 
1 5 band diagonal matrix, the matrix elements situated outside a predetermined band around the 
main diagonal of the approximated noise correlation matrix have the value zero, and 

the magnetic resonance image is reconstructed from the regularly resampled 
magnetic resonance signals which have been resampled from acquired magnetic resonance 
signals on the basis of the sensitivity profiles and the approximated noise correlation matrix. 
20 The magnetic resonance image is derived from the magnetic resonance signals 

sub-sampled in the k-space while utilizing the sensitivity profiles. Sub-sampling means that 
the sampling in the k-space is coarser, i.e. with a resolution in the k-space which is coarser 
than sufficient for the field of view of the magnetic resonance image. In a magnetic 
resonance imaging method the smallest wavelength of brightness variations in the magnetic 
25 resonance image relates to the field of view. The smallest wavelength is notably proportional 
to the magnitude of the field of view and to the sampling density in the k-space. In the case of 
sub-sampling the sampling is coarser than sufficient for the desired size of the field-of-view 
of the magnetic resonance image. The signal values are encoded on the basis of their wave 
vectors in the k-space and on the basis of the sensitivity profiles. The magnetic resonance 
30 signals of the respective receiving antennas correspond to respective signal channels. The 
noise contribution to the signals in each of the signal channels is a linear combination of 
noise contributions from the relevant signal channel and from (in principle) all other signal 
channels. The receiving antennae are, for example receiving coils that are sensitive to the 
magnetic resonance signals. Preferably, surface coils are used as the receiving antennas. Such 
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surface coils are arranged on the body of the patient to be examined and pick up notably 
magnetic resonance signals which are generated in the body of the patient to be examined in 
positions situated near the surface coil. The noise correlation between the signal channels is 
represented by a noise correlation matrix. For realistic numbers of magnetic resonance 
5 signals for a magnetic resonance image of diagnostic quality, the decoding of the magnetic 
resonance signals in the k-space and on the basis of the sensitivity profiles into pixel values 
for individual pixel positions in an image matrix constitutes, if no steps are taken, a matrix 
inversion problem requiring a high calculation capacity and long calculation times. The noise 
correlation ma tter may be approximated by the unity matrix, a block diagonal matrix 

10 or a two diagonal ma trix, which are all special examples of block diagonal or band 

diagonal matrices. Reconstruction of the magnetic resonance image from the sub-sampled 
magnetic resonance signals on the basis of the SENSE-algorithm includes optimisation of 
noise properties in the magnetic resonance image. This optimisation involves a noise 
correlation matrix which contains in the diagonal elements noise in the sampled magnetic 

1 5 resonance signals and in the off-diagonal elements noise correlations between respective 
sampled magnetic resonance signals acquired by different receiver antennae. It appears that 
as an approximation the noise correlation matrix may be replaced by the unity matrix. An 
alternative, more subtle, approximation is based on the recognition that the noise correlations 
are about constant over time. Hence, it appears that the noise correlation can adequately be 

20 described by a matrix having a sparse structure, i.e. approximately (block) diagonal. This 
sparse structure allows a virtual re-sampling or re-binning of the sub-sampled magnetic 
resonance signals from the respective receiver coils into virtual channels as linear 
combinations of the sub-sampled magnetic resonance signals from the individual receiver 
coils. The weights involved in this linear combination are obtained from the so-called 

25 Cholesky decomposition of the noise correlation matrix into a matrix product of an invertible 
left triangular matrix and its Hermitian conjugate. Then, the effective noise correlation matrix 
connecting the virtual channels is the unit matrix. It has been found according to the 
invention that in practice the correlation between noise contributions in the magnetic 
resonance signals from individual receiving antennas can be suitably approximated by means 

30 of a simpler matrix with contributions from only the vicinity of the main diagonal. It has even 
been found that this noise correlation may be replaced by the unity matrix. It has even been 
found that such a simplification strongly mitigates the matrix inversion problem, so that only 
a comparatively short calculation time and a limited calculation capacity are required. The 
magnetic resonance image can thus be reconstructed from the magnetic resonance signals 
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within a short period of time. It has also been found in practice that a 128x128 image matrix 
(soN= 128) can be reconstructed from the sub- sampled magnetic resonance signals within 
one minute. It has also been found that the approximation of the noise correlation matrix in 
the reconstruction of the magnetic resonance image on the basis of the sensitivity profiles 
5 does not have significant adverse effects on the diagnostic quality of the magnetic resonance 
image. This means that the magnetic resonance image has a suitable contrast resolution so 
that details of low contrast are suitably visibly reproduced in the magnetic resonance image. 
The reconstruction time of the magnetic resonance image is further reduced notably by 
applying a known fast Fourier transformation to the regularly resampled magnetic resonance 

10 signals. Regular (re)sampling means that the magnetic resonance signals in the k-space are 
sampled on a regular square grid. It has been found that for an N x N image matrix this 
simplification leads to a reduction of the matrix inversion problem from the order of 
magnitude N 4 to the order of magnitude N 2 or N 2 logN. 

The invention offers a high degree of freedom in choosing the acquisition 

15 trajectory to be followed through the k-space for the acquisition of the magnetic resonance 
signals. Such an acquisition trajectory according to the invention may give rise to irregular 
sampling of the k-space. It is notably not necessary to perform the sampling on a regular 
square grid in the k-space during the acquisition of the magnetic resonance signals. Thus, for 
example, respective parts of the k-space can be traversed at different speeds. The invention 

20 notably offers the possibility of choosing an essentially spiral-shaped trajectory through the 
k-space. Magnetic resonance signals are then acquired first from a central part of the k-space, 
so with wave vectors of comparatively small magnitude, after which magnetic resonance 
signals are acquired with a continuously faster increasing magnitude of the wave vector. Such 
an acquisition along a spiral-shaped trajectory, or a trajectory comprising one or more spiral- 

25 shaped steps in the k-space, is particularly suitable for use in MR angiography. Therein, 
magnetic resonance images of a patient to be examined are formed immediately after the 
administration of a contrast medium to the patient, for example an intravenously injected 
medium. The magnetic resonance signals from the center of the k-space relate mainly to 
rather coarse structures in the magnetic resonance image, including the arterial part of the 

30 vascular system of the patient to be examined. The venous part of the vascular system 
involves mainly much finer structures. When a spiral-shaped trajectory is followed, the 
magnetic resonance signals can be acquired from the arterial part before the contrast medium 
reaches the veins. Because, moreover, sub-sampling is applied, the acquisition of the 
magnetic resonance signals will not require much time. The combination of sub-sampled 
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acquisition of the magnetic resonance signals and scanning the k-space along a spiral-shaped 
trajectory enables fast acquisition of magnetic resonance images of the arterial part of the 
vascular system with a high spatial resolution. 

These and other aspects of the invention will be elaborated on the basis of the 
5 following embodiments which are defined in the dependent Claims. 

Preferably, respective receiving coil images are reconstructed from the 
magnetic resonance signals from individual signal channels, so from individual receiving 
antennas. Receiving coils are preferably used as the receiving antennas. Due to the sub- 
sampling of the magnetic resonance signals from the individual signal channels, aliasing 

1 0 artefacts such as fold-back phenomena occur in such receiving coil images.The 

reconstruction of the receiving coil images utilizes the approximated noise correlation matrix 
according to the invention. The magnetic resonance image is derived from the receiving coil 
images on the basis of the sensitivity profiles. The reconstruction of the magnetic resonance 
image on the basis of the receiving coil image and the sensitivity profiles is known per se as 

1 5 the SENSE method. This SENSE method is known per se from the article by Priissmann et 
al. in Proceedings ISMRM (1998), 579and from the article by Priissmann and Weiger in 
MRM42 (1999), pp. 952-962. The SENSE method enables sub-sampled acquisition of the 
magnetic resonance signals, thus reducing the time required for the acquisition of the 
magnetic resonance signals. 

20 It is alternatively possible to form magnetic resonance signals with full 

sampling from the sub-sampled magnetic resonance signals, while utilizing the sensitivity 
profiles, by combination of the sub-sampled magnetic resonance signals. The magnetic 
resonance image is reconstructed from the magnetic resonance signals obtained by 
combination. Various magnetic resonance signals in the k-space are then combined so as to 

25 fill in the lines in the k-space which have been skipped during the acquisition. This method is 
known by the acronym SMASH and is known per se from United States patent No. 
5,910,728. 

When receiving coils or surface coils are used as the receiving antennas, the 
coil sensitivity profiles of the receiving coils correspond to the sensitivity profiles of the 
30 receiving antennas. 

Preferably, the receiving coils are preferably decoupled essentially inductively. 
As the degree of inductive coupling of the receiving coils is lower, the noise level and noise 
correlation are lower. The noise level of the magnetic resonance image is thus reduced. 
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Preferably, an iterative inversion algorithm is used for the reconstruction of 
the magnetic resonance image. That is, the magnetic resonance image is reconstructed by 
iteration from the sub-sampled magnetic resonance signals. Starting with some initial vector, 
iterative algorithms yield a progression of approximate solutions, which converges to the 
exact solution. A variety of such techniques exist for the treatment of large linear systems. 
The so-called conjugate-gradient (eg) method is particularly suited. On the one hand, , it may 
be combined with FFT for very efficient calculations. On the other hand, the CG iteration 
does not require particular provisions for ensuring convergence. It converges safely given 
that the matrix involved is positive definite, which holds true for the matrix that connects the 
pixel-values of the reconstructed magnetic resonance image to the sub-sampled magnetic 
resonance signals through the gradient encodings and the coils sensitivity profiles. The CG 
algorithm theoretically yields the exact solution of an N 2 xN 2 system after at most N 2 
iterations. For N in the range of 128, though, it is not practical to carry out the entire 
procedure until mathematically strict convergence is achieved. However, in practice it 
appears that approximations leading to a good diagnostic quality for the reconstructed 
magnetic resonance image are obtained after a relatively small number of iterations already. 
Each CG iteration step consists in multiplying the matrix to be inverted with a residuum 
vector and several further calculations of minor complexity. Thus the iteration speed depends 
crucially on how fast matrix-vector multiplication can be performed. The number of 
iterations necessary to achieve a given accuracy is related to the so-called condition of the 
matrix to be inverted and the suitability of the starting vector. Because the dimensionality and 
size of the matrix inversion for the method according to the invention, iterative inversion 
algorithms are faster than, for example, direct inversion algorithms. Particularly favorable 
results are achieved, for example, by means of a Jacobi procedure, a Gauss-Seidel procedure 
or conjugated gradient (CG) methods. 

The invention also relates to a magnetic resonance imaging system suitable for 
carrying out the magnetic resonance imaging method according to the invention. A magnetic 
resonance imaging system according to the invention is defined in the independent Claim 8. 

The magnetic resonance imaging system according to the invention includes a 
control unit with a computer provided with a (micro)processor whereby the temporary 
gradient fields and RF excitations are controlled. The functions of a magnetic resonance 
imaging system according to the invention are carried out preferably by means of a suitably 
programmed computer or (micro)processor or a special purpose processor provided with 
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integrated electronic or optoelectronic circuits especially designed for the execution of one or 
more of the magnetic resonance imaging methods according to the invention. 

The invention also relates to a computer program with instructions for 
executing a magnetic resonance imaging method. It is a further object of the invention to 
provide a computer program whereby one or more of the magnetic resonance imaging 
methods according to the invention can be carried out. A computer program according to the 
invention is defined in the independent Claim 9. When such a computer program according 
to the invention is loaded into the computer of a magnetic resonance imaging system, the 
magnetic resonance imaging system can execute one or more magnetic resonance imaging 
methods according to the invention. Thus, on the basis of the instructions of the computer 
program according to the invention the technical effects can be realized for producing 
magnetic resonance images according to the invention. For example, a magnetic resonance 
imaging system according to the invention is a magnetic resonance imaging system whose 
computer is loaded with a computer program according to the invention. Such a computer 
program can be stored on a carrier such as a CD-ROM. The computer program is then loaded 
into the computer by reading the computer program from the carrier, for example by means 
of a CD-ROM player, and by storing it in the memory of the computer of the magnetic 
resonance imaging system. It is to be noted, however, that it is also possible to load the 
computer program according to the invention into the memory of the computer of the 
magnetic resonance imaging system via a network, for example the worldwide web. 

These and other aspects of the invention are apparent from and will be 
elucidated, by way of non-limitative example, with reference to the embodiments described 
hereinafter and the accompanying drawing. 

Brief designation of the Figures 

The Figure shows diagrammatically a magnetic resonance imaging system in 
which the invention is used. 

iii. Description of Figure 

The Figure shows diagrammatically a magnetic resonance imaging system in 
which the invention is used. The magnetic resonance imaging system includes a set of main 
coils 10 whereby the steady, uniform magnetic field is generated. The main coils are 
constructed, for example in such a manner that they enclose a tunnel-shaped examination 
space. The patient to be examined is slid into this tunnel-shaped examination space. The 
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magnetic resonance imaging system also includes a number of gradient coils 11,12 whereby 
magnetic fields exhibiting spatial variations, notably in the form of temporary gradients in 
individual directions, are generated so as to be superposed on the uniform magnetic field. The 
gradient coils 11, 1 2 are connected to a controllable power supply unit 2 1 . the gradient coils 
5 11,12 are energized by application of an electric current by means of the power supply unit 
21. The strength, direction and duration of the gradients are controlled by control of the 
power supply unit. The magnetic resonance imaging system also includes transmission and 
receiving coils 13, 16 for generating the RF excitation pulses and for picking up the magnetic 
resonance signals, respectively. The transmission coil 1 3 is preferably constructed as a body 

10 coil 13 whereby (a part of) the object to be examined can be enclosed. The body coil is 

usually arranged in the magnetic resonance imaging system in such a manner that the patient 
30 to be examined is enclosed by the body coil 13 when he or she is arranged in the magnetic 
resonance imaging system. The body coil 13 acts as a transmission antenna for the 
transmission of the RF excitation pulses and RF refocusing pulses. Preferably, the body coil 

15 13 involves a spatially uniform intensity distribution of the transmitted RF pulses (RFS). The 
same coil or antenna is usually used alternately as the transmission coil and the receiving 
coil. Furthermore, the transmission and receiving coil is usually shaped as a coil, but other 
geometries where the transmission and receiving coil acts as a transmission and receiving 
antenna for RF electromagnetic signals are also feasible. The transmission and receiving coil 

20 13 is connected to an electronic transmission and receiving circuit 15. 

It is to be noted that it is alternatively possible to use separate receiving coils 
16. For example, surface coils 16 can be used as receiving coils. Such surface coils have a 
high sensitivity in a comparatively small volume. The transmission coils, such as the surface 
coils, are connected to a demodulator 24 and the received magnetic resonance signals (MS) 

25 are demodulated by means of the demodulator 24. The demodulated magnetic resonance 
signals (DMS) are applied to a reconstruction unit. The receiving coil is connected to a 
preamplifier 23. The preamplifier 23 amplifies the RF resonance signal (MS) received by the 
receiving coil 16 and the amplified RF resonance signal is applied to a demodulator 24. The 
demodulator 24 demodulates the amplified RF resonance signal. The demodulated resonance 

30 signal contains the actual information concerning the local spin densities in the part of the 

object to be imaged. Furthermore, the transmission and receiving circuit 15 is connected to a 
modulator 22. The modulator 22 and the transmission and receiving circuit 15 activate the 
transmission coil 13 so as to transmit the RF excitation and refocusing pulses. The 
reconstruction unit derives one or more image signals from the demodulated magnetic 
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resonance signals (DMS), which image signals represent the image information of the 
imaged part of the object to be examined. The reconstruction unit 25 in practice is 
constructed preferably as a digital image processing unit 25 which is programmed so as to 
derive from the demodulated magnetic resonance signals the image signals which represent 
5 the image information of the part of the object to be imaged. The signal on the output of the 
reconstruction unit is applied to a monitor 26, so that the monitor can display the magnetic 
resonance image. It is alternatively possible for the magnetic resonance signals to represent a 
three-dimensional density distribution. Such a three-dimensional density distribution can be 
displayed on the monitor 26 in various ways; for example, projections to be selected by the 

10 user are displayed or stereoscopic pairs of images. It is alternatively possible to store the 
signal from the reconstruction unit 25 in a buffer unit 27 while awaiting further processing. 

The magnetic resonance imaging system according to the invention is also 
provided with a control unit 20, for example in the form of a computer which includes a 
(micro)processor. The control unit 20 controls the execution of the RF excitations and the 

15 application of the temporary gradient fields. To this end, the computer program according to 
the invention is loaded, for example, into the control unit 20 and the reconstruction unit 25. 

From the magnetic resonance imaging method which is known as the SENSE 
method the following relation is known to exist between the measured magnetic resonance 
signals (m) and the brightness or contrast values (/): 

20 

I = (E Hx ¥' ] Ey l E Hx ¥- l m (1) 

The encoding matrix E therein relates to the spatial encoding in the k-space on the basis of 
the wave vector of the magnetic resonance signals as well as to the spatial encoding by the 
sensitivity profiles (s Y )(r p ), where the index y represents the relevant surface coil and r p the 
25 position of the relevant volume element or voxel in the object to be examined. E denotes the 
complex adjoint of the encoding matrix. The noise correlation matrix is denoted as T, having 
the matrix elements 

r 

where a T represents the standard deviation of the noise in the signal channel x and co yT 
30 represents the weighing factor of the contribution to the noise in the signal channel y by the 
signal channel x. 

The encoding matrix E thus has the matrix elements: 
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Eyro ^ S y(. r n) Q KP 



(2) 



y,tc,p ~y\ -p 

Therein, s^r p ) denotes the 

spatial sensitivity profile at the area r p of the receiving antenna, notably the surface coil y. 
According to the invention the noise correlation matrix can be adequately approximated by 
5 a block diagonal matrix, or even by the unity matrix; this is notably so when the surface coils 
are inductively decoupled. The effective encoding is then reduced to 



I = {E H E)~ x E H m 



(3) 



10 In the more subtle approach it is assumed that noise correlation between different coils is 
constant over time as the underlying mechanisms are invariant. Receiver noise is then 
described by a sparse, i.e. block-diagonal matrix with a simple structure: 

^y,k;/y',k'; =VI/ Yy 5 K,K'- 

The time-independent, matrix *F can be determined experimentally by statistical analysis of 
15 reference noise samples taken in the absence of MR signal. Let r| y denote the noise output of 

the -th channel. Then the entries of *F are given by 



where the bar indicates time averaging. In this approach noise correlation of simultaneously 
acquired magnetic resonance signals is fully accounted for 
20 With such simplified noise statistics, the noise variance matrix can be eliminated by a simple 
trick. The basic idea is to create a set of virtual receiver channels by linear combination of the 
original ones, such that the virtual channels exhibit unit noise levels and no mutual noise 
correlation. Suitable weighting coefficients for this purpose are given by the inverse of the 
matrix L obtained by the Cholesky decomposition 

25 *¥ = LL H . 

Virtual sampling data with decorrelated unit noise are thus obtained from the original 
samples by 

m y%° rr ~ L ~ J h y' m Y k • The net coil sensitivities associated with the virtual channels 

y' 

are given likewise by 

30 ^rj-^^rj, 

v' 
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f f \ 



\ 



dk 



=ZWK^') Je-^-v) -k) 

r V V k J J 

leading to the modified encoding matrix 

E(y™Xp ~ e lk * rp Sy 6Corr (r p ) . With sample values and sensitivities modified in this 

fashion, the newly combined channels can be treated exactly like physical ones. The noise 
variance matrix of the virtual channels is equal to identity and can thus be omitted when 
Eq.[6] is reformulated for image reconstruction from the modified data. Note that by the 
transition to virtual channels the solution of Eq.[6] is not altered (see Appendix B). In 
particular, the optimisation of SNR is preserved. 

Both options, decorrelation on the one hand, and neglecting noise correlation on the other 
hand, lead to the same simplified formula. Dropping the superscript for decorrelation, and 

indicating again the approximated encoding matrix as E it reads in both cases (cf. eq. (3) 
above): 

{E H E) I = E H m . 



The approximated encoding matrix is now represented by the reference E . It has been 

found that this matrix inversion problem can be quickly numerically solved by iterative 
inversion. The data m are the signal values (amplitude and phase) of the magnetic resonance 
signals. The code in (3) includes an inverse Fourier transformation, so that 



/ = Q3 ] m 

Therein, Q is a generic representation of the residual codes, so apart from the 
inverse Fourier transformation. This can be more exactly written as vector and matrix 



(4) 



(E y)^ = \e->"«' 



r 



Zj p V r P^ r p-^ 

V P 



dr 



J 



f 



(£ w m) p =£i;(r,) Je 



-iknr. 



f 



V K 



dk 

J 



elements: 
and 
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For fast execution of the inverse Fourier transformations, a so-called "gridding" 
transformation is performed on the measured data on a regular square grid in the k- space: 

5 m = G®m, where G is, for example, a Gaussian convolution kernel: 



Therein, the inverse Fourier transformation can be quickly executed by means of known fast 
10 Fourier transformation (FFT) algorithms. The integrals in these expressions are efficiently 
calculated by FFT with gridding preparation using Kaiser-Bessel windows and two-fold 
oversampling. The computational cost of evaluating each of these equations is on the order of 



1 5 storing of E E . For the repeated evaluation of E Ey , the gridding approach can be 
taken advantage of in two ways: 

One can calculate and store the integral initially and perform common matrix-vector 

multiplication, creating the entries of E E row by row without storing them. This procedure 
20 has the advantage that it is compatible with the Gauss-Seidel-method, which often converges 
faster than the other methods. On the other hand the complexity of one matrix-vector 
multiplication remains essentially N 4 . 

IT * 

One can perform the evaluation of E E y in two successive gridding and FFT steps, This 

has the advantage that the cost for the matrix -vector multiplication is only on the order of N 2 
25 log N and that the procedure is easily implementable on parallel processing hardware. 

Furthermore, the magnetic resonance image is reconstructed from the residual code in the 
code Q by iterative inversion. Generally speaking, such iterative inversions can be written as: 

/(<+i) j0) + a /±( m 0) ^ ) 5 where the signal data m 0) is calculated on the basis of the 
30 residual code Q. from the reconstructed image data flK Each time a new estimate of the image 
data f l ^ J) is formed on the basis of the difference function A, which new estimate more 




(5) 



Therefore: 



G 



(6) 



O A 

only N log N as opposed to N for normal matrix-vector multiplication. The storing of the 

O A 

function f in eq. [5] requires memory size on the order of N instead of N for straightforward 
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closely fits each time the signal data measured by means of the magnetic resonance signals. 
The iteration commences with an estimate of the image data f 0J and the measured magnetic 
resonance signals m (0 K 
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1 . A magnetic resonance imaging method for forming a magnetic resonance 
image, wherein 

magnetic resonance signals are acquired by means of receiving antennas via a 
plurality of signal channels, 

which individual receiving antennas have respective sensitivity profiles, and 

the noise correlation between individual signal channels is represented by a 
noise correlation matrix, where 

the magnetic resonance signals are acquired with sub-sampling, 

regularly resampled magnetic resonance signals are resampled on a regular 
sampling grid from the acquired magnetic resonance signals, 

the noise correlation matrix is approximated by a block diagonal matrix or a 
band diagonal matrix, the matrix elements situated outside a predetermined band around the 
main diagonal of the approximated noise correlation matrix have the value zero, and 

the magnetic resonance image is reconstructed from the regularly resampled 
magnetic resonance signals which have been sampled from acquired magnetic resonance 
signals on the basis of the sensitivity profiles and the approximated noise correlation matrix. 

2. A magnetic resonance imaging method as claimed in Claim 1, wherein the 
approximated noise correlation matrix is a diagonal matrix, notably the unity matrix. 

3. A magnetic resonance imaging method as claimed in Claim 1, wherein 

respective receiving coil images are reconstructed from the regularly 

resampled magnetic resonance signals from individual signal channels, and 

the magnetic resonance image is derived from the receiving coil images and 
the sensitivity profiles. 

4. A magnetic resonance imaging method as claimed in Claim 1 or 2, wherein 
the sensitivity profiles of individual receiving antennas are essentially decoupled, from one 
another. 
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5. A magnetic resonance imaging method as claimed in one of the preceding 
Claims, wherein 

the magnetic resonance image is reconstructed from the regularly resampled 
5 magnetic resonance signals, notably from the receiving coil images, by means of an iterative 
inversion algorithm. 

6. A magnetic resonance imaging method as claimed in one of the preceding 
Claims, wherein 

1 0 - the magnetic resonance signals are acquired on the basis of a trajectory 

through the k-space which corresponds to a sampling of the acquired magnetic resonance 
signals outside a regular grid of sampling points in the k-space. 

7. A magnetic resonance imaging method as claimed in Claim 6, wherein the 
1 5 trajectory through the k-space includes an essentially spiral-shaped segment. 

8. A magnetic resonance system for forming a magnetic resonance image and 
arranged to 

acquire magnetic resonance signals by means of a plurality of receiving 
20 antennas and via a plurality of signal channels, 

which individual receiving antennas have respective sensitivity profiles, and 
the noise correlation between individual signal channels is represented by a 
noise correlation matrix, wherein 

the magnetic resonance signals are acquired with sub-sampling, 
25 - regularly resampled magnetic resonance signals are resampled on a regular 

sampling grid from the acquired magnetic resonance signals, 

the noise correlation matrix is approximated by a block diagonal matrix or 
band diagonal matrix, the matrix elements situated outside a predetermined band around the 
main diagonal of the approximated noise correlation matrix have the value zero, and 
30 the magnetic resonance image is reconstructed from the regularly resampled 

magnetic resonance signals which have been sampled from acquired magnetic resonance 
signals on the basis of the sensitivity profiles and the approximated noise correlation matrix. 



9. 



A computer program which includes instructions for 
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resampling regularly resampled magnetic resonance signals on a regular 
sampling grid from the acquired magnetic resonance signals, wherein 

the acquired magnetic resonance signals are acquired with sub-sampling via a 
plurality of signal channels, and 

the noise correlation between the individual signal channels is represented by a 
noise correlation matrix, 

the noise correlation matrix is approximated by a block diagonal matrix or a 
band diagonal matrix, the matrix elements situated outside a predetermined band around the 
main diagonal of the approximated noise correlation matrix have the value zero, and 

the magnetic resonance image s reconstructed from the regularly resampled 
magnetic resonance signals which have been sampled from acquired magnetic resonance 
signals on the basis of the sensitivity profiles and the approximated noise correlation matrix. 
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